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ABSTRACT 



We study the process pp — > J/tp + 7 + X at the Tevatron energy (y/s = 1.8 TeV). 
The perturbative QCD contributions to this process from the gluon-fusion and the 
fragmentation mechanisms are computed. For the entire range of p? that can be 
probed at the Tevatron, the fusion contribution is found to be dominant. Consequently 
the QCD prediction for this process has better precision than the purely hadronic 
production of J/ip. An experimental study of this process at the Tevatron will help to 
clarify several aspects of quarkonium production at large pr- 
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The production of quarkonia at large pt is usually described by the parton-fusion 
mechanism in the framework of the colour-singlet model |], [|, f|. In this model, the 
wave-function at the origin for a S'-wave (or its derivative for a P-wave) quarkonium 
is convoluted with the cross-section for producing a heavy quark pair with the proper 
spin, parity and charge-conjugation assignments. The wave-function |Po(0)| 2 can be 
estimated from the S'-wave quarkonium decay widths within an uncertainty of ~ 30%, 
arising from relativistic correction 0, ||. The corresponding uncertainty in the esti- 
mate of the derivative of the wave- function IP'^O)! 2 from P-wave quarkonium decay 
is much larger, since the NLO correction for this case is not yet available [[§, fj|, f|. 
This uncertainty propagates into the hadroproduction cross-section of J /if), which is 
dominated by the x ~ > V* + 7 decay contributions. Together with the uncertainty from 
the choice of QCD scale in the cc pair production, one can predict the J /if) hadropro- 
duction cross-section to within a factor of ~ 3. Within this normalisation uncertainty, 
the colour singlet model prediction agrees with the J /if) hadroproduction data from 
fixed-target and ISR experiments over a reasonable range of p? 0, [|, 0, • A compre- 
hensive account of the model prediction and comparison with these data is given in 



At higher energies where 6-quark production becomes important, the decay of b 
quarks is a mechanism that contributes to J /if) production, in addition to the parton- 
fusion mechanism of the colour- singlet model ||. Recently, the CDF experiment stud- 
ied the large-pr J /if) production at the Tevatron collider and, using the secondary 
vertex information, separated out the 6-quark decay contribution. The remaining J/^'s 
are presumably produced directly by the fusion mechanism, but the CDF experiment 
found that the cross-section is an order of magnitude larger than that predicted by 
the fusion mechanism of the colour- singlet model. Even if one takes into account the 
above mentioned normalisation uncertainty in the theoretical predictions, it does not 
suffice to restore agreement with the CDF data. 

The anomalously large cross-section measured in the CDF experiment - the CDF 
effect - indicates that there is a new mechanism through which J /if) production takes 
place at large pr- At high-energy colliders, where gluons and charm quarks are co- 
piously produced, the fragmentation of these particles into J /if) has been suggested 
as an important mechanism for large-pr J/ip production pD[ . The validity of this 
suggestion is borne out by the explicit computation of the fusion and fragmentation 
contributions for the Tevatron energy, presented in Refs. [jTTJ, 12, [H|. In these papers, 



it was shown that the fragmentation contribution is the dominant contribution at the 
Tevatron energy; and together with the fusion contribution it can explain the large 
cross-section measured by the CDF experiment to within a factor of 2-3. The largest 
contribution comes from gluons fragmenting into %'s, which subsequently decay into 
J/^'s. In Ref. [lT|, the energy dependence of the fragmentation and fusion contribu- 



tions was studied as well. It was shown that at ISR energies the fusion contribution 



dominates over the pt range of interest, although the inclusion of the fragmentation 
contribution helps to improve the quantitative agreement with the ISR data. 

In Refs. [|11], 113, the production of ip' by fragmentation was also studied. Unlike 
J/ip production, ip' production does not involve the P-state contribution. Conse- 
quently, the fragmentation contribution is comparatively smaller. The total predicted 
ip' cross-section is atleast an order of magnitude smaller than the cross-section mea- 
sured by CDF. This large discrepancy with ip' data seems to indicate a large contribu- 
tion to the S'-wave quarkonium production from a still unknown mechanism. As noted 



in [11], [12| , this will show up in a larger enhancement of the ip' cross-section than the 
J/ip. Whatever be the source of this new contribution, it is important to probe it phe- 
nomenologically through other processes of quarkonium production at large pr- In this 
note, we point out that the production of a large-pr J/4> with an associated photon 
will serve as a very interesting probe for such mechanisms of quarkonium production. 

As we shall see below, the associated production of a large-p^ J/ip with a photon 
has three distinctive features in comparison with the purely hadronic J/ip production 
discussed above. 1) The fusion contribution to this process dominates over the frag- 
mentation for the entire pt range of interest at the Tevatron energy. 2) The fusion 
contribution comes entirely from the direct production of the 5-state. 3) The S-state 
wave function |i?o(0)| 2 is known to reasonable precision; and besides the uncertainty 
associated with the QCD scale is much smaller here than the purely hadronic process. 
In short, the dominant contribution to this process comes from the direct produc- 
tion of the S-state, which can be predicted to a reasonable precision. Thus it is well 
suited to probe for an anomalous enhancement of the .S-wave quarkonium production 
at large-pr- 

In recent years, there have been several discussions of associated production of a 



J/ip and a photon in hadron-hadron experiments |JI1, IR [TO, O . This production of a 



J/ip—^f pair can also take place at HERA [[14, l8|, through resolved photon processes. In 



the framework of the colour- singlet model, the production of a J/ip and an associated 
photon, proceeds through the following subprocess: 

g + g^J/iP + 7 . (1) 

C-invariance forbids a x coupling to two gluons and a photon; so a \ cannot be 
produced via the above subprocess. Of course, the x contribution through gg — > 
X — > J/4> + 7 is possible; but this will contribute to \ow-pr J/^'s and 7's, and can 
be eliminated by a pt or invariant mass cut. Therefore, large-pr J/ip+'J production 
through gluon fusion involves only the S"-wave resonance. 

In addition to the gluon-fusion contribution, one expects a fragmentation con- 
tribution to the associated production of J/ip + 7 at large-p^- The fragmentation 



contribution to J/tp is expected to come mainly from the gluon and charm-quark jets 



10| , Pi [TJ, |13[- Thus the basic process of interest is the associated production of a 
large-pT photon with a gluon (or a charm quark), followed by the fragmention of the 
latter into J/tp either directly or via the x' s - Consequently the fragmentation contri- 
bution involves both S and P states. Another distinctive feature of the fragmentation 
contribution is the emergence of the J/tp along with the other jet fragments. Thus 
in principle one could distinguish it from the fusion process (1) by requiring either 
J/tp isolation or p? balancing between J/tp and 7. It would be hard to implement 
this in practice, however, since the J/tp is expected to carry away the bulk (~ 70%) 



of the jet momentum [12|. Therefore one has to deal with the sum of the fusion and 



fragmentation contributions. 

There is another source of J/tp + 7 production at large-pr - i-e. the associated 
production of a large-p^ J/tp with a hadron jet, followed by the fragmentation of the 
latter into a hard photon [[B5], ^, [2T| . However this process can be substantially reduced 



by the standard isolation cut on 7 used in the direct photon experiments. It was noted 
in [^T|] that even with the isolation cut one should not neglect this fragmentation 
contribution while dealing with the NLO correction to direct photon production. In 
the present work we shall work with the lowest-order direct photon production cross- 
section and hence neglect the contribution from photons produced via fragmentation. 

Finally, we wish to make an observation here concerning the model-dependence 
of the fusion contribution. We will use the colour-singlet model of J/tp production in 
this work; but there exists another model of J/tp production - the semi-local duality 
model [22]. In this model, the J/tp cross-section is given by integrating the open cc 



cross-section between 2m c and the open cc production threshold and multiplying by 
a normalisation factor, which is not specified by the model but has to be determined 
by comparing with data. In this model, the cc pair produced in the hard-scattering is 
not required to be a colour-singlet, neither is a projection to the quarkonium quantum 
numbers required. The colour and quantum numbers are assumed to rearrange them- 
selves appropriately, in the process of soft-gluon emission. Consequently, in this model, 
the gg-fusion production of J/tp + 7 is allowed, in addition to the gluon-gluon fusion 



16| . However, at the large energies of our present interest the gg fusion is expected to 
dominate in the duality model, so that the predictions of the two models would not be 
significantly different. 

The subprocess cross-section for J/tp+j production in the colour- singlet model is 
given as 



da, T/ , . 167raa^M|i? (0)| : 



dt y ™ /r " 27s 2 
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where s, i and u are the usual Mandelstam variables and M is the mass of the J/ip. 
To obtain the hadronic pt distribution, this subprocess cross-section is folded in with 
the gluon densities, as follows: 

<la -{AB^J/^ + 1 ) = 



dp T 

Y, dy dx 1 xig A ( y x 1 )x2g B (x 2 )- - — --^(99 -»• J/i/> + l)- (3) 

J J 2x\ — xre y at 

In the above equation, g& and g B are the gluon distributions in A and B, x\ and x 2 

are the respective momentum fractions. Energy-momentum conservation determines 

x 2 to be 

xix T e~ y — 2r 

X2 = 2 Xl -x T ev ' (4) 

where r = M 2 /s, y^s the centre-of-mass energy, and y the rapidity at which the J/ip 
is produced. We also have 

2M T 2p T 



x T = Jx T + 4r = — ^, x T = —;=■ (5) 



As discussed above, in addition to the fusion process there is a contribution to 
J/ip+'-f production from the fragmentation of gluons and charm quarks. At lowest 
order, the subprocesses that give rise to a parton at large pt with an associated photon 
are the usual direct photon subprocesses: 

q + q -»• 9 + 1 
Q + 9 ^Q + l, (6) 

where we have used Q to denote the charm quark, and the charm-quark density in the 
initial state is that generated by the Q 2 evolution of the structure functions. Assuming 
factorisation, the fragmentation contribution to the J/ip+'-f cross-section can be written 

as 

da(AB -+ (J/V>, Xi) + 7 + X) = J2 j dz da(AB -> c + 7 + X)D c ^ (J/ ^ Xi) (z, fi), (7) 

where da(AB — > c + 7 + X) is the hard-scattering cross-section for the production of a 
parton (denoted by c) in association with a photon, D(z,fi) is the fragmentation func- 
tion specifying the fragmentation of parton (gluon or charm quark) into the required 
charmonium state, and z, as usual, is the fraction of the momentum of the parent 



parton carried by the charmonium state. The sum in the above equation runs over all 
contributing partons. The fragmentation function is computed perturbatively at an 
initial scale /io, which is of the order of m c . It is then evolved to the scale typical of the 
fragmenting parton, which is of the order of Pt/z, using the Altarelli-Parisi equation. 
The full set of initial fragmentation functions needed to obtain the J/ip and the \ 
contributions have now been computed. They are D g ^j/^ |10[, D g ^ x p3j, P c ^v fH 



and D c ^ x f25| , 



For the fragmentation process, the cross-section is given by a formula similar to 
Eq. |H but with an extra integration over z, or equivalently over x 2 - We have 

^-(AB^(J/iP, Xi ) + 1 X) = 
dp T 

Y^ J dydx 1 dx 2 G a /A(xi)G b /B(x2)D c ^j/ i)ai )(z) — --ir( a b -> c + 7), (8) 
with z given by 

2 y x 2 x x ! y ! 

For da/dt(ab — > c + 7), we have used the lowest-order direct photon cross-sections [20 . 



In Fig. 1, we have presented our results for the fusion and fragmentation contri- 
butions to J/ip + 7 production in pp collisions at y^ =1.8 TeV. We have plotted 
Bda/dpx as a function of px, integrated over a pseudo-rapidity range \i]\ < 0.5, where 
B is the J/ip branching ratio into leptons (B = 0.0594). In our computations, we have 
used pTj the updated MRSD-' parametrisations |28j for the parton densities in the 
nucleon. The parton densities are evolved to a scale Q 2 = /i 2 /4, where /1 is chosen to 
be Mt for the fusion case, and equal to p^ c = Pt/z for the fragmentation case. The 
fragmentation functions are evolved to the scale Pt/z. Another uncertainty that en- 
ters the normalisation of the predicted cross-sections is that due to the wave-function 
at the origin, R , the derivative for the P-states, R[ (or equivalently Hi) and H' 8 , 
which is a parameter that describes the g — > x fragmentation through a colour-octet 
mechanism. The last parameter is related to the infrared divergent term in the NLO 
correction to the P-state decay and has by far the largest uncertainty. For the param- 
eters Po, Hi and Pg, we have used the values quoted in Refs. [Tt], |2"3| : Pq = 0.8 GeV 3 , 
Hi = 15.0 MeV, H' 8 = 3.0 MeV. 

The fusion contribution (shown as the solid line in Fig. 1) comes entirely from 
the Estate production and hence depends only on Pg. Being a gluon-gluon fusion 
process, it is seen to dominate over the whole range of px considered. The gluon 
fragmentation contribution (including g — > J/ip and g — > x an d shown as the dashed 
line in Fig. 1) and the charm-quark fragmentation (shown as the dotted line) are more 
than an order of magnitude smaller than the fusion contribution, even at the largest 



values of pr considered. These features may be contrasted with the purely hadronic 
J/ip production where the gluon fragmentation dominates over fusion and both the 



contributions are dominated by the production of P-states [11, 12]. 



In Fig. 2, we show the J/ip+'j cross-section (i.e. the sum of the fusion and frag- 
mentation contributions) as a function of p?- This is shown as the solid line in Fig. 2. 
The cross-section is large enough to be measurable at the Tevatron. We also estimated 
the magnitude of the ip'+j cross-section, but it turns out to be two to three orders of 
magnitude smaller than the J/ip+j cross-section, and may be too small to be observed. 
We have also considered the effect of varying the scale from /i/2 to 1\x on the J/ip + 7 
cross-section (the latter choice is shown as the dashed line in Fig. 2). The cross-section 
is remarkably stable under scale variation. This is due to a substantial compensation 
between the decrease of a s and the increase of the low-x gluon distribution g(x ~ .01) 
in eqs. 2,3 as the scale is increased from /x/2 to 2/x. As a result the QCD prediction 
for the dominant process of gluon fusion (1) is remarkably stable under scale variation 
at the Tevatron energy. The net uncertainty from the choice of scale as well as the 
QCD parametrisation may be roughly estimated as ±40%. As remarked before, the 
uncertainty from the estimate of the corresponding wave function R^ is also quite small 
(~ ±30%). Combining the two gives an overall uncertainty of ±50% for the predicted 
J/tp ± 7 cross-section of Fig. 2. Measurement of the cross-section at the Tevatron will 
therefore be very useful in probing for any new source of quarkonium production at 
large px- In particular an anomalous enhancement of the S'-wave quarkonium produc- 
tion would show up here as clearly as in the purely hadronic production of ip' if it is 
associated with the fusion mechanism, but not if it comes from the fragmentation. 

In summary, we have analysed the associated production of a large-p^ J/ip and 
7 at the Tevatron energy. The perturbative QCD predictions for the fusion and the 
fragmentation contributions have been computed. The fusion contribution is seen 
to dominate this cross-section over the entire pt range of interest, unlike the purely 
hadronic production of J/ip. More over the fusion contribution to this process comes 
entirely from the direct production of the S'-wave quarkonium. Consequently one has 
a rather precise QCD prediction for this process. Experimental measurement of this 
cross-section will help to clarify several aspects of quarkonium production at large pr, 
and in particular those related to the anomalous hadronic J/ip and ip' cross-sections 
recently reported by the CDF experiment (the CDF effect). 

We would like to thank Drs. P. Sphicas and R.S. Bhalerao for discussions. 
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Figure captions 

Fig. 1 Lower figure: The cross-section Bda/dpT (integrated over the pseudorapidity 
range —0.5 < r] < 0.5) for the process pp — > J/ip + 7 + X as a function of Pt 
at y^i =1.8 TeV. The different curves correspond to the direct production via 
fusion (solid line), the gluon fragmentation contribution (dashed line) and the 
charm-quark fragmentation term (dotted line). 

Fig. 2 The scale dependence of the J/ifi + 7 cross-section as a function of px for pp 
collisions at y^ =1.8 TeV. The solid curve in both figures is for the scale /i/2 
and the dashed curve is for the scale 2/x. 



